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Abstract.  A  spherical  condensed  phase  (volatile  particle)  is  placed  in  an  infinite  expanse  of  a  binary  mixture  consisting  of 
the  vapor  of  the  condensed  phase  and  another  noncondensable  gas.  The  mixture  is  assumed  to  be  at  rest  with  uniform  pressure 
and  temperature  at  infinity.  However,  the  vapor  and  the  noncondensable  gas  have  small  uniform  density  (or  concentration) 
gradients  of  the  same  magnitude  but  in  the  opposite  direction.  The  flows  of  both  components  around  the  particle  are  analyzed 
numerically  by  a  finite-difference  method  on  the  basis  of  a  model  Boltzmann  equation  for  gas  mixtures,  and  the  force  acting 
on  the  particle  (diffusiophoretic  force)  is  obtained  accurately  for  a  wide  range  of  the  Knudsen  number. 


INTRODUCTION 

Thermophoresis,  photophoresis,  and  diffusiophoresis  are  typical  phenomena  due  to  the  effect  of  gas  rarefaction  that 
manifest  themselves  in  connection  with  aerosol  particles  and  have  been  important  subjects  in  kinetic  theory  of 
gases  [1,  2,  3].  Recent  advances  in  microtechnology  have  reinforced  their  importance.  Among  these  phenomena, 
diffusiophoresis  [4,  5,  6,  7,  8],  the  force  acting  on  an  aerosol  particle  and  its  resultant  motion  in  a  mutually  diffusing 
background  mixture  of  gases,  is  peculiar  to  gas  mixtures  and  thus  has  more  nonintuitive  features. 

In  the  present  study,  we  investigate  the  diffusiophoresis  of  a  volatile  particle.  To  be  more  specific,  we  consider 
a  spherical  condensed  phase  placed  in  a  binary  mixture  at  rest  consisting  of  the  vapor  of  the  condensed  phase  and 
another  noncondensable  gas,  each  of  which  has  a  small  density  (or  concentration)  gradient  of  the  same  magnitude  but 
in  the  opposite  direction.  In  this  situation,  evaporation  and  condensation  of  the  vapor  taking  place  on  the  surface  of 
the  particle  may  have  an  important  effect  on  the  flows  of  both  components  around  the  particle.  We  will  analyze  such 
flows  for  a  wide  range  of  the  Knudsen  number  numerically  on  the  basis  of  the  linearized  version  of  a  model  Boltzmann 
equation  [9,  10]  and  obtain  the  force  acting  on  the  particle  (diffusiophoretic  force)  accurately. 


FORMULATION  OF  THE  PROBLEM 
Problem  and  assumptions 

Let  us  consider  an  infinite  expanse  of  a  binary  mixture  of  two  gases,  gas  A  and  gas  B.  There  is  no  macroscopic  motion 
of  the  total  mixture,  and  its  temperature  and  pressure  (and  thus  its  molecular  number  density)  are  uniform.  But.  the 
molecular  number  densities  of  respective  gases  have  uniform  gradients  with  the  same  magnitude  but  in  the  opposite 
directions.  Let  us  denote  the  uniform  temperature  by  To,  the  uniform  pressure  by  po,  and  the  uniform  molecular 
number  density  by  ?iq  (=  po  j kTq,  where  K  is  the  Boltzmann  constant),  and  let  us  express  the  number  density  of  gas 
A  as  +  (dnA /dX\  j^A)  and  that  of  gas  B  as  Uq  —  (dnA /dX\  )oo2fi ,  where  =  no,  and  Xj  is  a  rectangular  space 

coordinate  system.  Correspondingly,  the  concentration  of  gas  A  is  expressed  as  Xq  +  (<?/t  /dX\  fWi ,  and  that  of  gas 
B  as  Xq  ~  {dXA  /  dXfj^Xi,  where  Xq  =  «o/«o,  Xq  =  «f/«o  (thus,  Xq+Xo  =  1),  and  {dxA/dXi)aa  =  nfl  {dnA / dXf)oo. 

In  such  a  mixture  of  gases  (say,  at  the  origin  A)  =  0),  a  spherical  particle  of  radius  L  and  temperature  To,  made  of 
the  condensed  phase  of  gas  A,  is  placed.  Therefore,  evaporation  or  condensation  of  gas  A  (vapor)  may  take  place  on 
the  surface  of  the  particle.  We  investigate  the  steady  behavior  of  the  vapor  (gas  A)  and  the  noncondensable  gas  (gas  B) 
around  the  particle  and  obtain  the  force  acting  on  the  particle  (diffusiophoretic  force)  for  a  wide  range  of  the  Knudsen 
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number  on  the  basis  of  kinetic  theory  under  the  following  assumptions: 

(i)  The  behavior  of  the  mixture  is  described  by  the  model  Boltzmann  equation  proposed  by  Hamel  [9]  or  that  by 
Garzo  et  al.  [10].  [Once  the  linearization  around  an  equilibrium  state  at  rest  is  made,  these  two  models  give  the 
same  equation.  We  will  actually  handle  the  linearized  equation,  as  mentioned  in  (v)  below.] 

(ii)  The  molecules  of  the  vapor  leaving  the  surface  of  the  particle  are  distributed  according  to  the  Maxwellian 
distribution  corresponding  to  the  vapor  in  the  saturated  equilibrium  state  at  rest  with  temperature  Tq. 

(iii)  The  molecules  of  the  noncondensable  gas  undergo  diffuse  reflection  on  the  particle  surface. 

(iv)  The  is  equal  to  the  molecular  number  density  of  the  vapor  in  the  saturated  equilibrium  state  mentioned  in  (ii). 

(v)  The  imposed  gradient  of  the  number  density  (or  the  concentration)  of  each  gas  is  small  [\{dnr / dX\)„\  -C  n^/L, 
or  \  (dxA / dX\)oo  \  <IC  Xo  /L],  so  that  the  equation  and  boundary  condition  can  be  linearized. 

The  assumption  (iv)  as  well  as  the  assumption  that  the  temperature  of  the  particle  is  the  same  as  that  of  the  mixture 

at  infinity  (To)  is  not  essential  to  the  present  result,  as  will  be  explained  in  “Results  and  discussions”.  The  case  where 

there  is  a  uniform  flow  of  the  mixture  at  infinity  will  also  be  discussed  there. 


Basic  equations 

We  first  summarize  the  main  notations.  In  what  follows,  the  Greek  letters  a  and  /3  stand  for  the  labels  A  and  B 
of  the  component  gases.  For  gas  a  (a  =  A,  B),  Fa  is  the  velocity  distribution  function,  na  the  molecular  number 
density,  v“  the  flow  velocity,  Ta  the  temperature,  and  ma  the  mass  of  a  molecule.  Our  reference  state  is  the 
equilibrium  state  at  rest  with  temperature  To  and  number  densities  of  respective  gases  and  nB ,  which  is  expressed 
as  F0“  =  (2kTq / mA)~3 /2 Ea  with  Ea  =  (/?z“/7r)3/2exp(— maC,f),  where  ma  =  ma  /mA  and  is  the  dimensionless 
molecular  velocity  [(2KTo/mA)l'2Q  is  the  dimensional  one].  Then,  we  let  Fa  =  F0“(l  +(f>a),  na  =  «“( 1  +  <»“), 
v“  =  (2fc7o/MiA)1/2Mi“,  Ta  =  7o(l  +  t“),  so  that  \(pa\,  |<ua|,  |m“|,  |t“|  <C  1.  In  addition,  x,  =  X,/L  and  nff  =  «“  j nAy 
For  the  model  equation,  the  collision  frequency  of  the  a -gas  molecules  for  the  collisions  with  the  /3-gas  molecules  is 
given  by  Ka^n^  with  constants  Ka&  (KAB  =  KBA).  Then,  the  mean  free  path  of  the  vapor  molecules  (/  in  the  reference 
equilibrium  state  is  given  by  =  (8KTo/7tinA)1/2(KAAny  +  KAB  nB)~l .  With  KaP  and  £q,  we  define  KaP  =  KaP / K AA 
and  k  =  (yTr/ 2)Kn=  (y/n/2 )(£q/L),  where  Kn  is  the  Knudsen  number. 

With  the  above  notations,  the  linearized  model  equation  is  written  as 


d<t>L 


1 


£  -  (/)“),  ( a=A,B ), 


^  dXi  (1  +KABnB)kptA.l 

y^ccp  =(0a+  2 Aa  ^.uaP  +  {^aQ  _  3/2j  %ap , 

hif+mPu?  ^ ap  =za+  2mamP 
ma  +  mP  ’  (ma  +  mP)2 


a(3 

U;  = 


(t*3  -t“), 


J  <f>aEad 3C,  uf  =  J  Q(j)aEad2t;,  Za  =  (2/3)  J (maQ  -  3/2 )<j)ccEccd2’ £ , 


(1) 

(2) 

(3) 

(4) 


where  d 3^  =  dC\d ^d ^3 .  The  boundary  condition  on  the  surface  of  the  particle  (|jc,-|  =  1)  is  given  as 


</>A  =  0,  </>b  =  -2 (nmB)1/2  f  £jnj$BEBd3t;,  (£■«,•  >0),  (5) 

0 

where  77,  is  the  unit  normal  vector  to  the  particle  surface  pointing  into  the  gas.  On  the  other  hand,  the  boundary 
condition  at  infinity  (|x,  |  — >  <x>)  is  given  by 


</>A  — ►  </>A  =  Coo(xi  —  2^i  Dk), 


=  -^(xi-2^Dk), 
no 


~  (l+kABnB)(\+mB) 
2KAB(\+mBnB) 


(6) 


where  c<»  =  (L/% o)(dxA / dX\)oo  =  (L/n^)(dnA / dX\)oo,  and  </)“  is  the  solution  to  Eq.  (1)  expressing  the  situation 
described  in  the  first  paragraph  in  “Problem  and  assumptions”.  That  is,  this  solution  gives  the  following  macroscopic 
variables: 


coA  =  Co 


,xi,  coB  = —CooXi/nB,  t“  =  0,  7/f  =  —  CooDk,  uB  =  CooDk/mBnB ,  =  w  “  =  0. 


(7) 
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FIGURE  1.  Profiles  of  the  macroscopic  quantities  of  the  vapor  for  various  values  of  k  in  the  case  where  mA  =  mB ,  =  KAB  = 

Kbb ,  and  Hq  =  nB.  (a)  ft>A,  (b)  r4,  (c)  ,  (d)  zfg. 


Note  that  there  is  a  flow  of  each  gas  (i.e.,  diffusion)  in  the  direction  opposite  to  its  density  or  concentration  gradient 
though  there  is  no  flow  of  the  total  mixture  [ u \  =  (uA  +  mB  nBuB)  /  (\  +  mBhB)  =  0,  where  (2kt7o/wa)1/2m!-  is  the  flow 
velocity  of  the  total  mixture] . 

It  is  seen  from  Eqs.  ( 1) — (6)  that  we  need  to  specify  the  following  five  dimensionless  parameters:  mB ,  KAB,  KBB ,  nB, 
k  to  solve  the  problem  (we  do  not  need  to  specify  Coo  because  of  the  linearity  of  the  problem). 

NUMERICAL  ANALYSIS 

In  the  case  of  a  single-component  gas,  spherical  geometry  allows  a  convenient  similarity  solution  to  the  linearized 
Boltzmann  equation  [11,  12].  In  the  present  problem,  the  similarity  solution  of  the  same  form,  i.e., 

4>“  =  C,  0?)cos0  +  C o*?(r,  C,  e?)sin0,  (8) 

is  shown  to  be  compatible  with  Eqs.  (1)— (6),  where  (r.  6,  <p)  is  the  spherical  coordinate  system  with  its  origin  at  the 
center  of  the  particle  and  its  polar  axis  along  the  x\  axis,  Cr,  Co,  and  C,p  are,  respectively,  the  r,  0,  and  <p  components 
of  Ci,  C  =  (C/2)1/2’  and  %  is  the  angle  between  Ci  and  the  r  direction  [  A  -  =  cos  1  ('Qr !  'Q)^  0  <  6'  <  n\.  Thus,  the 
independent  variables  are  reduced  to  r,  £,  and  6U .  Correspondingly,  coa/ cos  6,  w“/cos0,  »g/sin0,  and  rr//  cos  0 
turn  out  to  be  functions  of  r  only,  where  w“,  Ug,  and  u’^  (=  0)  are  the  r,  9,  and  <p  components  of  it" . 

We  solve  the  boundary  value  problem  for  and  <f>“  numerically  by  a  finite-difference  method.  The  method 
is  basically  the  same  as  that  used  for  a  single-component  gas  in  Ref.  [12]  where  a  slow  rarefied  gas  flow  past  a 
sphere  is  analyzed,  except  that  we  are  handling  a  simpler  model  collision  term  here  rather  than  the  original  Boltzmann 
collision  integral.  The  difficulties  that  we  still  share  with  Ref.  [12]  are:  (a)  the  solution  approaches  the  state  at  infinity 
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(a) 


CQ 


(b) 


r 


very  slowly;  (b)  the  velocity  distribution  functions  exhibit  discontinuities  around  the  sphere.  These  are  overcome  by 
essentially  the  same  way  as  in  Ref.  [12].  That  is,  (a)  is  solved  by  connecting  the  numerical  solution  with  the  explicit 
asymptotic  solution  for  the  far  field.  The  discontinuities  in  (b)  are  described  accurately  by  a  combination  of  a  standard 
finite-difference  method  and  a  method  of  characteristics  that  was  devised  in  Ref.  [13]. 

As  a  reference,  we  give  the  data  about  the  computational  systems  that  are  used.  For  mB  =  0.5,  the  computational 
domain  is  1  <  r  <  vm  \xm  =  200 k  ( k  <  5),  750  (k  =  5),  and  1000  (k  =  10)],  0  <  £  <  8,  and  0  <  6^  <  n,  and  the  numbers 
of  the  grid  points  in  r,  £,  and  dr  are  601  (nonuniform),  71  (nonuniform),  and  401  (uniform),  respectively.  The  systems 
for  other  mB  are  the  same  except  that  0  <  £  <  7  and  61  grid  points  (nonuniform)  are  used  for  C- 


RESULTS  AND  DISCUSSIONS 

First  we  show  the  profiles  of  the  macroscopic  quantities.  Figure  1  shows  the  profiles  for  the  vapor,  that  is, 
coA/ Coo  cos  8  —  r ,  TA/cooCos0,  uA  / c^  Dk  cos  6.  and  Ug/cooDk  sin  6  (note  that  uf  — >  —c.„Dk  as  r  — >  °°)  versus  r,  for 
various  k  in  the  case  where  nr4  =  mB,  KAA  =  KAB  =  KBB ,  and  =  nB  (or  mB  =  1,  KAB  =  KBB  =  1,  and  hB  =  1). 
Figure  2  shows  the  corresponding  profiles  of  the  noncondensable  gas,  that  is,  n/’/rv.  cos  6  +  r/nB,  xB /crx  cos  6, 
mBnBuB /cooDk  cos  6,  and  mB nB uB / CooDk  sin  8  (note  that  uB  — >  CooDk/mBnB  as  r  — >  °o)  versus  r,  in  the  same  case. 
Let  us  suppose  that  cTO  >  0.  Then,  there  is  a  uniform  flow  of  gas  B  in  the  x\  direction  at  infinity.  Therefore,  the  flow 
velocity  field  of  gas  B  around  the  particle  is  more  or  less  similar  to  that  in  the  case  of  a  uniform  (noncondensable)  gas 
flow  past  a  sphere  [12],  On  the  other  hand,  there  is  a  uniform  flow  of  gas  A  in  the  negative  x\  direction  at  infinity,  so 
that  the  flow  velocity  field  of  gas  A  has,  to  some  extent,  a  common  feature  with  that  in  the  case  of  a  uniform  vapor 
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FIGURE  3.  The  hd  versus  k  in  the  case  where  KM  =  KAB  =  KBB .  (a)  mB /mA  =  0.5.  (b)  mB /mA  =  1,  (c)  mB /rtA  =  2.  The  result 
for  the  free-molecular  flow  is  shown  by  a  solid  line  (;ig/«Q  =  0.1),  a  dashed  lone  (//q/«q  =  1),  and  a  dot-dashed  line  (n£/nfi  =  10). 


flow  past  a  spherical  condensed  phase  [14].  Thus,  condensation  takes  place  on  the  head  (the  part  facing  to  the  positive 
x\  direction)  of  the  particle  and  evaporation  on  the  tail  (the  part  facing  to  the  negative  x\  direction)  [see  Fig.  1(c)].  The 
evaporation  and  condensation  are  enhanced  by  the  imposed  density  gradient  c0 „  that  makes  the  number  density  around 
the  head  higher  than  the  saturation  number  density  of  the  vapor  (at  the  particle  temperature)  and  that  around  the  tail 
lower  than  the  saturation  number  density.  The  flow  field  around  the  sphere  is  roughly  a  superposition  of  these  effects. 

Next,  we  show  the  result  of  the  force  acting  on  the  particle  (diffusiophoretic  force).  The  force  Ff  =  ( Ff ,  0,  0)  can 
be  expressed  in  the  following  form: 


Ff 


1/2 

PqL~Dab 


hd, 


_  (mA  +  mR)  kTq 
AB~  mAmBKABn0  ’ 


(9) 


where  Dab  is  the  mutual  diffusion  coefficient  for  the  present  model  Boltzmann  equation.  The  hd  versus  k  is  shown 
in  Fig.  3  in  the  case  where  KAA  =  KAB  =  KBB  =  1.  Figure  3(a)  is  for  mB /mA  =  0.5,  Fig.  3(b)  for  mB /mA  =  1,  and 
Fig.  3(c)  for  mB /nf4  =  2.  In  each  figure,  the  results  for  nB/n q  =0.1,  1,  and  10  are  shown.  The  hd  is  positive  in  most 
cases.  However,  it  is  rather  sensitive  to  the  molecular  mass  ratio  mB /mA  and  becomes  negative  when  mB /mA  =  2  and 
k  is  relatively  large.  The  result  for  the  free-molecular  flow  (k  =  °°): 

hd  =  (4^/3)Vm®5[8(l  -  Vot®)  +  7t]/(1+7hs),  (10) 

where  D  is  given  in  Eq.  (6),  is  also  shown  in  Fig.  3. 

In  the  present  study,  we  assumed  that  the  temperature  of  the  particle,  say  Tw,  was  the  same  as  that  of  the  gas  at 
infinity  (To)  and  that  the  reference  number  density  of  the  vapor  (;/q)  was  equal  to  the  saturation  number  density,  say 
/if,  of  the  vapor  at  temperature  Tw  (=  To).  If  these  conditions  are  not  satisfied  (Tw  is  still  assumed  to  be  uniform),  but  if 
the  deviations  are  small  enough  for  the  linearization  to  be  valid  (|7o  —  Tw\/Tq  <S  1  and  \nf  —  nA  \ /  nf  <<C  1),  then  we  only 
need  to  add  a  spherically  symmetric  term  of  the  form  <f>“ ( r,  C.  6^)  to  the  right-hand  side  of  Eq.  (8),  which  does  not 
contribute  to  the  force  acting  on  the  particle.  Therefore,  Eq.  (9)  is  also  valid  in  this  more  general  situation.  Incidentally, 
a  uniform  Tw  is  realized  when  the  thermal  conductivity  of  the  particle  is  high. 

So  far,  we  have  considered  the  case  where  there  is  no  flow  of  the  total  mixture  at  infinity.  When  there  is  a  slow 
uniform  flow  Uc„j  =  (I/oci,  0,  0)  of  the  total  mixture  at  infinity  [|t/ooi  |  «;  (2fc7b/'«A)1/2],  we  need  to  add  2£i«ooi  and 
2mB £iMoci  to  (l>A  and  (f>B  in  Eq.  (6),  respectively,  where  uc„ i  =  (2kTq / mA )_ 1  2f/._oi .  For  instance,  if  we  consider  the 
case  where  the  noncondensable  gas  is  at  rest,  we  let  u^_  \  =  —CooDk/mBnB.  Because  of  the  linearity  of  the  problem, 
the  solution  of  such  a  problem  can  be  obtained  as  the  sum  of  the  solution  presented  in  this  paper  and  the  solution 
of  the  problem  of  a  uniform  flow  of  the  mixture  past  a  spherical  condensed  phase,  that  is,  the  problem  in  which  the 
same  spherical  condensed  phase  is  placed  in  the  uniform  equilibrium  flow  of  the  mixture  with  temperature  7o,  velocity 
Uooi ,  and  number  densities  of  respective  gases  //()  and  nB.  We  have  also  analyzed  this  problem  by  the  same  numerical 
method.  The  drag  force  Ff  =  {Ff ,  0,  0)  acting  on  the  particle  in  this  case  is  expressed  in  the  following  form. 

Ff  =  bnpLU^hD,  jU  =  KTo^/ir^K^  +n%KAB)  +n%/ {r^KBA +  n%KBB)\,  (11) 
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FIGURE  4.  The  ho  versus  k  in  the  case  where  =  KAB  =  KBB .  (a)  mB /mA  =  0.5,  (b)  mB /mA  =  1,  (c)  mB /mA  =  2. 


where  u  is  the  viscosity  coefficient  for  the  present  model  equation.  The  hD  versus  k  is  shown  in  Fig.  4  for  KAA  = 
Kab  =  KBB .  In  summary,  the  force  acting  on  the  particle  when  there  is  a  uniform  flow  of  the  total  mixture,  in  addition 
to  the  density  gradient  of  each  component  considered  originally  in  this  paper,  is  given  by  Ff  +  Ff. 

Finally,  let  us  go  back  to  the  original  problem  without  a  flow  of  the  total  mixture  at  infinity.  If  the  particle  is  left 
without  constraint,  it  starts  moving  because  of  the  force  (9)  (diffusiophoresis).  The  final  velocity  of  the  particle  can  be 
obtained  by  using  the  result  in  the  preceding  paragraph.  That  is,  we  superpose  at  infinity  a  uniform  flow  (say  [/,■)  of  the 
mixture  that  vanishes  the  total  force  (Ff  +  Ff  =  0).  Then  the  velocity  of  the  particle  relative  to  the  flow  of  the  total 
mixture  (—[/,)  gives  the  final  velocity  of  the  particle.  From  Eqs.  (9)  and  (11),  we  obtain  the  velocity  of  the  particle 
vf  =  {Vf,  0,  0)  as 


yd  p()L  /  2k~7q  \  1/2  Dab  (  dxf\  hd 

1  6n  y  mA  J  fi  \  dXi  )  x  ho 


(12) 


Following  Ref.  [15],  we  can  show  that  the  unsteadiness  caused  by  the  particle  motion  through  the  mixture  with  a 
nonuniform  concentration  can  be  neglected  because  of  the  smallness  of  the  concentration  gradient. 
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